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Abstract 

NMR spectroscopy is a excellent tool for monitoring in-situ chemical reactions. 
In particular, DOSY measurement is well suited to characterize transient species 
by the determination of their sizes. However, here we bring to light a difficulty 
in the DOSY experiments performed in out-of-equilibrium systems. On such a 
system, the evolution of the concentration of species interferes with the mea¬ 
surement process, and creates a bias on the diffusion coefficient determination 
that may lead to erroneous interpretations. 

We show that a random permutation of the series of gradient strengths used 
during the DOSY experiment allows to average out this bias. This approach, 
that we name p-DOSY does not require changes in the the pulse sequences nor 
in the processing software, and restores completely the full accuracy of the mea¬ 
sure. This technique is demonstrated on the monitoring of the anomerization 
reaction of a- to /3-glucose. 

Keywords: DOSY; out-of-equilibrium; complex mixtures; permutation; 
glucose 


Introduction 

Diffusion Ordered SpectroscopY (DOSY) is an NMR technique separat¬ 
ing components in a mixture spectroscopically according to their translational 
diffusion [HE]. Considered as an NMR chromatographic technique, DOSY NMR 
has been widely used in life sciences and chemistry for the analysis of complex 


* Corresponding author 

Email address: madelsucOunistra.fr (Marc-Andre Delsuc) 
^ These two authors contributed equally 





mixtures, such as biological extracts [3], euviroumeutal fluids [4], supramolecu- 
lar chemistry [5l [ 6 ], uauoparticles[3 [ 8 ], or pharmaceutical formulatious[9]. For 
spherical shaped compoueuts the diffusiou coefhcieut D is iuversely proportioual 
to the hydrodyuamic radius vh as described by the Stokes-Eiusteiu equatiou: 


D = 


ksT 

67rr]rH 


( 1 ) 


where is the Boltzmauu’s coustaut, T is the absolute temperature, aud 77 
is the viscosity. Other shapes cau be haudled either with exteuded aualytical 
models [To] or by a statistical descriptiou pTl [T^ . 

A DOSY experimeut cousists iu a series of spiu echo or stimulated echo spec¬ 
tra recorded with iucreasiug gradieut field streugth. Wheu measured with the 
Stimulated Echo pulse sequeuce (STE), diffusiou coefficieuts cau be determiued 
by aualyziug the sigual amplitude as a fuuctiou of the square of the gradieut 
pulse area usiug the Stejskal-Tauuer equatiou [13] giveu here for squared pulses 

I = (^2) 

where Iq is the sigual iuteusity obtaiued iu the abseuce of field gradieut aud / 
iu the preseuce of a gradieut of streugth D is the diffusiou coefhcieut, 6 is the 
gradieut pulse duratiou, A is the ehective dihusiou delay dehued as the time 
betweeu the start of the bipolar gradieut pairs, aud 7 is the gyromaguetic ratio 
of the spiu uuder study; see [14] for details 

The aualysis of the dihusiou coefficieuts is a powerful tool to gaiu iusight 
ou molecular sizes aud shapes [HI \TE\ [TB] . aud cau also be used to iuvestigate 
molecular iuteractious [T7HT9] . polymer polydispersity|20[[2T], or to characterize 
reactive iutermediates [22H25] . 

However, wheu measuriug DOSY experimeut ou au evolviug system, a dif- 
hculty arises. Accordiug to the classical Stejskal-Tauuer equatiou, the iuteusity 
at each gradieut streugth depeuds ou the iuteusity Iq with uo gradieut applied. 
Usually, lo is measured at the begiuuiug of the DOSY experimeut aud cou- 
sidered coustaut throughout the measuremeut. However, iu out-of-equilibrium 
systems, Iq is uot coustaut auymore but varies coutiuuously, due to evolutiou 
of the couceutratious. Cousequeutly, throughout the DOSY measuremeut, the 
sigual iuteusity measured for a giveu species will depeud ou the couceutratiou 
evolutiou as well as the iucreasiug gradieut streugth, aud will appear to decay 
either too rapidly or too slowly, depeudiug ou the evolutiou of the couceutra¬ 
tiou, decreasiug or iucreasiug over the time of the measuremeut, creatiug a bias 
iu the aualysis. 

The kiuetic rates of the out-of-equilibrium system has to be takeu iuto ac- 
couut aud be compared to the timeframe of the DOSY measuremeut. Rapid 
DOSY measuremeuts mouitoriug systems with slow kiuetic rates are uot affected 
to a siguificaut degree, because ouly miuor couceutratiou chauges occur over the 
experimeutal time. Ou the other haud, loug DOSY measuremeuts mouitoriug 
systems with fast kiuetic rates are uot affected either, but fail to record dyuamic 
iuformatiou of the system such as reactive iutermediates. 
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The effect is most obvious for systems with intermediate kinetic rates mon¬ 
itored by DOSY measurements performed in a time similar to the kinetic con¬ 
stant of the chemical systems. In this case, the DOSY spectra display a shift of 
the peaks along the diffusion coefficient axis leading to wrong results. 

Nilsson et al.[26] proposed the use of trilinear multivariate analysis in order 
to characterize both the reaction kinetics and the diffusion coefficient, in the 
analysis of a set of diffusion experiments measured during the chemical reaction. 
In this case the concentration variations is used as a third dimension in the 
multiway analysis. However in this approach, the DOSY measurements are 
supposed to be rapid compared to the kinetic rates, and there is no compensation 
for the unavoidable bias in the diffusion coefficient measure. To our knowledge, 
the DOSY studies performed so far on out-of-equilibrium systems [22ll26] have 
not taken into account the analytical bias caused by concentration variations. 

We propose here to measure, within a single DOSY measurement, a series of 
stimulated echo spectra with pulse gradient strengths sampled in a permutated 
manner, as it was already proposed in order to reduce Fi artifacts in 2D NMR 
spectroscopy [27] . 

A permutation of the list of gradient strengths allows the separation of the 
experimental time, during which the system evolves, from the gradient strength 
sampling, thus removing the experimental bias. Knowing the actual evolution of 
the concentration, an optimal permutation could eventually be designed. How¬ 
ever this knowledge is usually not available, and such an optimal permutation 
not accessible, so we choose to perform random permutations of the gradient 
list. This permutation distributes randomly the concentration variation over 
the whole range of gradients, actually transforming an experimental bias into 
an additional random noise (see figure S14 in Supp. Mat.). The analysis of 
the permutated experiment leads to a less precise measure (larger error-bars) 
but eventually more accurate (less systematic error). Each realization of the 
random permutation may remove only partially the systematic error, and the 
accuracy is improved with a larger number of gradient levels. Simulations per¬ 
formed for several model systems (Figures S13 to S27 in Supp.Mat.) show that 
large bias may be observed for a regular approach and are fully reduced by this 
procedure. Additionally, the improvement of the measure with larger gradient 
lists is clearly observed. 

From the sequential list of pulse field gradient strengths, a random permu¬ 
tation is computed, and the DOSY is measured with a new list of gradient 
strengths in permutated order. This protocol has been named permutated- 
DOSY (p-DOSY). To test the effect of the permutated acquisition on the accu¬ 
racy of the measurement, a set of simulations was performed and the diffusion 
coefficient of glucose was monitored during the equilibration from pure a-glucose 
to the anomeric equilibrium of a and /3-glucose [28] (see Figure [^. 

Materials and Methods 

Simulations were performed in Python 2.7, using the libraries numpy 1.9 
and scipy 0.15. They are presented in supplementary material, and can be 
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Figure 1: Anomerization of a-Glucose to /3-glucose with proton numbering 


downloaded at https://github.com/delsuc/p-DOSY, 

Crystalized a-D-glucose (anhydrous, 96% purity) was purchased to Siguia- 
Aldrich (USA). Deuterated solvents DMSO-d6 (99.96%) and D 2 O (99.9%) were 
purchased to Euriso-Top (France). All samples were measured in 5.0mm tubes. 
A stock solution of 100.0 mM a-D-glucose in DMSO was prepared. 100 jjL of the 
glucose stock solution were added in 900 fiL D 2 O resulting in final concentration 
of 10.0 mM. 

DOSY NMR measurements were carried out both on a Bruker and on 
a Varian instrument, (for Varian experiments see supplementary material). 
The Bruker instrument was a 500 Avance I NMR spectrometer, operating at 
500.137 MHz for ^H, equipped with a 5 mm TXI probe. The pulse sequence used 
was a longitudinal eddy current delay bipolar gradient pulse (ledbpgp2s) An 
exponential gradient list of 32 values was created by using the standard AU pro¬ 
gram dosy. For the p-DOSY experiment, the list was permutated by an ad-hoc 
python script (see supplementary material), and the permutated list was used 
for acquisition and processing. Experiments were acquired with 8 scans, 5/2 oi 
3.0 ms and A of 150 ms, for an overall time of 23:09 min. A series of 15 spectra 
was acquired at a temperature of 300 K for a total duration of 6.0 hours. Before 
and after every DOSY spectrum, an 8 scans proton spectrum was recorded. 
The whole experiment was duplicated once with regularly increasing series of 
gradients, and once with a permutated series. 

The DOSY spectra acquired on Bruker were processed with the NMRnote- 
book program (version 2.7, NMRTEC France), using the DOSY MaxEnt imple- 
mentation[29j |30]. The experiments performed on Varian were analyzed with 
the vnmrj (version 6.3) software. In both cases, no modification had to be made 
to the computation protocol despite the permutated entries in the DOSY file. 

The final diffusion coefficients were calculated from the average position of 
the anomeric signal for each species. Error bars were estimated as the stan¬ 
dard deviation along the diffusion axis of the different signals. Minor viscosity 
variations were compensated by normalizing all measurements to the DMSO 
diffusion coefficient. 


Results and discussion 

The p-DOSY approach, and its capacity to restore accuracy in the deter¬ 
mination of diffusion coefficients was tested using glucose anomerization as a 
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Figure 2: Overlay of two standard DOSY experiments, acquired between: a) out-of- 
equilibrium (0.00-25.00 minutes) (blue color) and b) at equilibrium (5.41-6.05 hours) (light 
blue color) 


test case for out-of-equilibrium systems. Glucose crystalizes in the a form, as 
a consequence, this form is dominant in the solution just after dissolution. The 
anomerization of a to /3-glucose was followed by acquiring DOSY measurements 
for approximately 6.0 hours, with a-D-glucose varying from 96.6% to 44.6%. 
The same experiment was performed twice, the first time with the standard 
DOSY method, and the second time using the p-DOSY approach. 

The diffusion coefficients in water of both species a and ;d-glucose are ex¬ 
pected to be very similar, as they have very similar shapes. However, they 
appear with quite different diffusion coefficients when measured with a stan¬ 
dard DOSY experiment, acquired in a sequential manner, as can be seen in 
figures and 1^ In particular, the signals from the ;d-glucose are found initially 
at a markedly slower diffusion coefficient and recover the normal value over the 
course of the experiment. This is the result of the progressive increase of con¬ 
centration of /3-glucose over the duration of a single DOSY measurement, which 
creates a bias in the exponential analysis of the signal. At the end of the 6 hours 
experiment, the recovery is not complete. 

On the other hand, this bias is completely absent when the anomerization 
equilibrium is monitored using the p-DOSY measurements, run under the same 
conditions (see figure [^. 

In both experiments, the diffusion coefficients of water and DMSO display 
the standard values for all experiments run during the kinetics, indicating that 
the error on /3-glucose is not due to an experimental bias but solely due to the 
variation of the concentration. Figure presents the evolution of the DOSY 
signal intensities and of the apparent diffusion coefficients over the course of the 
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Figure 3: Overlay of two p-DOSY experiments, acquired between: a) out-of-equilibrium 
(0.00-25.00 minutes) (green color) and b) at equilibrium (5.37-6.00 hours) (yellow color) 


experiment. It shows that the initial /3-glucose concentration is very low, and 
that the two populations equilibrate slowly, so that 6 hours are not sufficient 
to reach the complete equilibrium. When measured with the standard DOSY 
experiment, the apparent ^d-glucose diffusion coefficient, varies in an anomalous 
manner before converging toward the correct value, while the p-DOSY measure 
presents stable measurements for both species. 

The poor accuracy of the diffusion coefficient of /3-glucose is attributed to 
miscalculation of the analysis of the reference intensity, J^. The diffusion coef¬ 
ficient of /3-glucose consequently, appears to have an incorrect low value. This 
effect is stronger in the first DOSY measurements since there is a higher kinetic 
rate. Similarly, the a-glucose that evolves in the opposite direction, displays 
an apparent diffusion coefficient slightly larger than expected in the first DOSY 
measurements. The effect is less marked because the relative concentration 
variation of the a-glucose is less important than in the case of /3-glucose. 

On the other hand, the experiment performed with a series of p-DOSY 
measurements does not display this systematic error, and accurate values are 
determined throughout the whole experiment. Because of the dispersion effect 
of the random permutation operation, the systematic evolution of the concen¬ 
tration is transformed here into random errors on the Iq value. This added 
random noise has no effect on the analysis, but disperses in a random man¬ 
ner the concentration variation over the gradient list, thus providing a way to 
average out bias. However, because of the added noise, a slightly less precise de¬ 
termination of the diffusion coefficients is expected. This effect, while not very 
strong, can be observed on Figure by the evolution of the error bars on the 
p-DOSY determined diffusion coefficients. The error bars present a somewhat 
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Figure 4: a) concentration evolution of ct, /3-glucose over experimental time, measured as 
the intensity of the anomeric proton peak based on ID experiments interleaved during the p- 
DOSY experiment series; b) evolution of the apparent diffusion coefficients of a and ^-glucose 
over experimental time based on a series of sequential conventional DOSY measurements; c) 
same as b, obtained with the p-DOSY experiment. In b) and c) the horizontal doted line is 
located at the diffusion coefficient of a-glucose, measured as the mean value of p-DOSY value 
over the whole period. 
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larger spreading in the beginning of the kinetic and converge toward a minimal 
value. 

In a nutshell, one can say that the p-DOSY experiment provides on evolving 
chemical systems, a slightly less precise, but much more accurate determination 
of the diffusion coefficient that the standard DOSY experiment. 


Conclusion 

The accurate estimation of diffusion coefficients is challenging in systems 
away from equilibrium. An inherent characteristic of out-of-equilibrium sys¬ 
tems is the evolution of the concentration of the different components during a 
DOSY NMR experiment. Since the basis of a single DOSY measurement is the 
analysis of the exponential decay of the signal with increasing gradient strength. 
Concentration evolution, within the characteristic measurement time, creates a 
bias, which can lead to large errors in the diffusion coefficient estimation. This 
was observed for the test case of the anomerization of glucose. Randomizing the 
list of gradient strength arrays within a single DOSY measurement, decouples 
the gradient evolution from the concentration evolution and removes the bias 
in the analysis, thus restoring the accuracy of the measure. The concentration 
evolution is still there though, as a noise source, and may impact the precision 
of the measure compared to a static sample, although its impact is minimal. 
Depending on the spectrometer software, this experiment requires little or no 
modification of the standard procedure, and the produced data-set can easily 
be analyzed using regular processing programs. 

Altogether, the DOSY experiment with a randomized list of gradients (the 
permutated DOSY or p-DOSY) is an unbiased experiment in the presence of 
an evolving system, in contrast to conventional DOSY. It is perfectly suited for 
the analysis of chemical reactions and out-of-equilibrium systems and can be 
used to monitor the rate of reaction kinetics, to characterize in-operando tran¬ 
sient species and measure their diffusion coefficients, or to observed molecular 
organization phenomena and dynamics effects. 

The robustness of p-DOSY can also be used to protect diffusion results from 
any experimental artifact, such as temperature shift or spectrometer drift when 
such artifacts are known and cannot be easily compensated for. As it does not 
add any burden in the acquisition step nor at the processing step, we recommend 
spectroscopists to use p-DOSY in this experimental cases. 


Supplementary Information 

Supplementary Information is available with this work. Figures S1-S2 present 
additional glucose spectra. Figures S3-S6 compare DOSY and p-DOSY mea¬ 
sured on a Varian spectrometer. Figure S7 presents the python program used 
for computing permutations. Figures S8-S27 present the results of the various 
simulations, the siumlation program can be downloaded at https: //github. 
com/delsuc/p-DOSY, 
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1. Glucose study by NMR 

The a,(3-D-glucose was studied by NMR and the assignment of the molecule 
have been done, to be able to characterize the a a and j8 protons for further 
analysis. This assignment is shown in the 1D spectrum in Figure S1 and it 
corresponds to the numbering of both anomeric molecules in Figure 1. 



Figure S1: Proton assignment of a mixture of a,p-D-giucose (42.7% a-D- 
glucose and 57.4%) of (S-D-giucose) at equiiibrium 
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2. p-DOSY 


The DOSY experiment and the p-DOSY experiment, measured on a-glucose 
sample after solubilisation at the beginning of the anomerization kinetics are 
shown superimposed in Figure S2. The bias on the diffusion coefficient can 
be seen for the 
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Figure S2: Overlay of two DOSY measurements, acquired between 0.0-25.0 
minutes, by p-DOSY arrays (green color) and sequential (blue color) DOSY 
arrays 


3. Varian experiments 

The Varian instrument was an Inova 500 NMR spectrometer operating at 
499.7722 MHz for ^H, equipped with a four-nucleus 5mm ^H(^®N-^^P) PFG 
High-Field Indirect detection probe. A DOSY bipolar pulse pair stimulated 
echo with convection compensation (Dbppste_cc) was used. 

Except for a temperature of 298.0 K, the experimental set-up equivalent to the 
Broker experiments. For the p-DOSY experiment, the permutated pulse 
gradient strength arrays was obtained by activating the randomization flag. 
DOSY experiments were run in four different modes: sequential, interleaved, 
permutated, and permutated and interleaved combined. 

The sequential DOSY spectra were acquired with 32 gradient 
increments with increasing gradient strengths. In the interleaved experiment, 
the block size was set to 4 (4 sets of 32 increments of 4 scans each). For the 
permutated experiment, 32 increments recorded by permutating the gradient 
strength pulses. A combination of both interleaved and permutated was also 
acquired (4 sets of 32 increments of 4 scans each with permutated gradient 
pulses). 
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Figure S3: Sequential DOSY measurement acquired between 0. 0-25.0 
minutes (13.0-50.0% ^-D-giucose) in Varian 



Figure S4: Interleaved DOSY measurement acquired between 0.0-25.0 
minutes (13.4-42.5 % p-D-glucose)in Varian 



Figure S5: P-DOSY measurement acquired between 0.0-25.0 minutes (8.3- 
40.8% (S-D-glucose) in Varian 



Figure S6: Combination of permutated and interleaved DOSY measurement 
between 0.0-2.0 (7.34-33.77%i ^-D-glucose) in Varian 
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4. Permutated Difframp and difflist 

Here is presented the "Permute" script to permutate the Difframp and difflist 
for acquisition and processing the p-DOSY experiments. 

The original Difframp, created by TopSpin when a DOSY experiment is 
launched, has to be replaced by the permutated Difframp before running the 
p-DOSY experiment. The same procedure has to been followed for the difflist 
before processing the spectrum. 


N = 64 # N is the length of the gradient list 

# Adapt to your experiment 
perm = random.sample (range (N), N) 

print N, "points perm 

filelist = ['difflist’, ’Difframp’] # difflist is for processing, Difframp is for acquisition 
for fname in filelist: 

F = open (fname, ’r’ ) 

FF = open(fname+’P’, ’w’) 

11 = [] 

for 1 in F: 

if not 1.startswithC ’#’ ): 

11.append(1) 
elif 1 != ’##END=’: 

FF.write(l) 

F. closeO 
if len(ll) != N: 

raise Exception( ’This prgm is meant for 7.d points ’7.N ) 
for i in range(N): 

FF.write(11[perm[i]]) 

FF.write( ’##END=’ ) 

FF.close 0 

python permute.py script used to compute a random permutation of the 
initial your Difframp and difflist files, for acquisition and processing of the 
p-DOSY experiment. 
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